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Heath-bath effects in the dynamics of atom + cavity system are studied. The temperature effects 
are explored using thermofield dynamics formalism. It is found that the dynamics of the system is 
sensitive to small changes in the temperature and the thermal effects lead to increasing instabilities 
by causing transitions from regular to chaotic motion. 



I. INTRODUCTION 



Cavity quantum electrodynamics deals with studying the interaction of atoms with photons in high-finesse 
cavities in a wide range of the electromagnetic spectrum, from microwaves to visible light. The fact that 
the system "atom + cavity mode" is a quantum system makes cavity quantum electrodynamics (QED) an 
excellent testing ground for such important issues of modern quantum physics as quantum measurement theory, 
entanglement, quantum computation, quantum interference and at the same time provides a unique possibility 
for trapping, cooling and manipulating of atoms 2].. Practical importance of cavity QED is mainly related 
to potential possibility for manipulating atoms and photons in mesoscopic scales. Therefore, in recent years 
cavity QED has become one of the hot topics both in theoretical and experimental context [l|43|- In particular, 
a number of new phenomena, such as the realization of quantum phase gate |4[, creation of Fock states of 
the radiation field 0, Q] and quantum non-demolition measurements Q have been found. The dynamics of 
a single atom trapped in a microcavity is governed by quantum electrodynamics. This makes cavity QED 
an interdisciplinary area as many subfield of physics, such as quantum and atomic optics, cold atom physics, 
physics of nanosized systems and quantum information, may use important results of the cavity QED. 

Recently cavity QED is considered in dealing with nonlinear dynamics 0, Q ■ Mapping quantum equations 
of motion onto classical ones, for the Jaynes-Cummings Hamiltonian, which includes recoil motion of the atom, 
Prants. et. al, explored phase-space dynamics of the atom interacting with a single cavity mode by analyzing 
Poincare surface sections and calculating Lyapunov exponents. In such a case classical and not quantum chaos 
described by classical equations of motion corresponding to Jaynes-Cummings hamiltonian are studied that are 
classical counterpart of the atom+photon+cavity system that are described by Prants et al This system 

can be considered as a periodically driven quantum system as most of the experiments on the study of chaos 
in atoms explore interaction of cold atoms with the standing waves of various type, such as phase-modulated 



standing wave 22j, amplitude-modulated standing wave and pulsed standing wave |23| 



In this paper we extend the approach developed by Prants et al. 0, Q to the case of the system coupled 
to a heat bath i.e. at finite temperature. For this purpose we use a real-time finite-temperature quantum 



field theory, thermofield dynamics |13Hl5j . In particular, by solving temperature-dependent classical equations 
of motion we explore Poincare surfaces sections, Lyapunov exponents and Levy flights in finite-temperature 
cavity QED. The role of a thermal bath, noise and dissipation effects in cavity QED is of importance both 
from fundamental and practical viewpoints. This paper is organized as follows: in the next section we present 
brief description of TFD formalism. Section 3 briefly recalls the results for T — obtained by Prants et al. 
Section 4 presents treatment of "atom + confining mode" system at finite temperature, while numerical results 
for Poincare surfaces of section, Levy flights and Lyapunov exponent at T — are given in section 5. Finally, 
section 6 presents some concluding remarks. 
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II. THERMOFIELD DYNAMICS 



TFD is a real time operator formalism of quantum field theory at finite temperature in which any physical 
system can be constructed from a temperature-dependent vacuum which is a pure state f]~3j] — fsoj] . For finite 
temperature it has been recognized for a long time that the Hilbert space has to be doubled. This is achieved 
with Thermofield Dynamics. One of the objective was to get a theory that real time and at the same time 
it is at finite temperature. The second Hilbert space is introduced along with a set of operators that are 
similar to be distinct from the normal operators. The final results don't depend explicitly on this second set of 
operators. However, Bogoliubov transformation mixes the two set of operators. This brings in factors that are 
temperature-dependent. The vacuum is also temperature-dependent and the creation and annihilation operators 
are converted to temperature-dependent operators, such that the annihilation operator acting on this vacuum 
gives zero. For a long time it was thought that these operators are ghosts. However, Umezawa showed that 
these operators can be viewed as degrees of freedom of the heat bath in the classical theory Later on 

it was shown by Santana and Khanna that one set gives the set of observables and the other gives symmetry 



operators. This argument is based on group theory and is given in the ref. 20 1 

Thus TFD is a powerful tool for exploring quantum dynamics of a system at finite temperature provided its 
Hamiltonian can be represented in terms of annihilation and creation operators. It has found many applications 



in condensed matter physics [13|, [16J, |17| , especially in superconductivity theory and related topics. Recently 
TFD has been applied to explore quantum chaos in the Yang-Mills-Higgs system [3] and for the calculation 
of the spectra of a strongly interacting bound system j^~oj ] . In this work we apply TFD prescription to Jaynes- 
Cummings model. It should be noted that TFD has been applied earlier to Jaynes-Cummings model (21] ] where 
the thermal noise effects in quantum optics are studied. 

Applying TFD prescription to a quantum system implies performing two actions: [bll . 

i) doubling of the Fock space which means that all operators arc doubled by introducing their tilded partners 
describing the effects of a heath bath; and 

ii) using Bogoliubov transformations, which makes Hamiltonian of the system temperature-dependent . With 
respect to a Hamiltonian operator written in terms of annihilation and creation operators, doubling means that 
the total Hamiltonian is written as 

H = H — H, (I) 
while the Bogoliubov transformations which are given by 

a = a((3) cosh 9 + a f (/3) sinh 9 

at = a f (/3) cosh + a(/3) sinh 9 

a = a t (/3)sinh0 + a(/3)cosh0 (2) 
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fit = a(J3)wi)iO + 5t(J3)coBh.e, 

make this Hamiltonian temperature-dependent. 
Here 

/3=-r^-; stib?6 = (ef'-l)- 1 (3) 
k B l 

and the annihilation and creation operators satisfy the following commutation relations: 

[oC8),o+C8)]=l [o(/9),o+( / 9)] = l (4) 
All other commutation relations are zero. 

We note that we are dealing with a system in equilibrium, i.e. the temperature, T, is constant. Since the 
Hamiltonian of the Jaynes-Cummings model is written in terms of annihilation and creation operators, it is 
convenient to use TFD formalism to treat this model at finite temperature. Applying the above prescription 
to the Jaynes-Cummings model we transform the Hamiltonian of the system into the temperature-dependent 
form that allows us to treat chaos using the same approach as that used by Prants et.al Q. 



III. CAVITY QED AT T = 

In this section we briefly recall the case of T = which is recently explored in detail in a series of papers 
by Prants et.al. Q-Ql- The simplest cavity QED system is a single two-level atom interacting with a single 
standing wave mode moving along the x-axis with the frequency U)f. 

The dynamics of this system in the presence of atomic recoil motion is described by Jaynes-Cummings 
Hamiltonian, which is written as 

P2 . 

H = 1- huj a S z + hajfa+a — hn {a+S- + aS + ) cos(kjx), (5) 

where S z , 5+ and <SL - are expressed in terms of pauli matrices i.e. = (S x ,S y ) and S z is the z - 

component, and a+ and a are the creation and annihilation operators, respectively, describing a selected mode 
of the radiation field of the frequency ujf and the wave number kf in a lossless cavity. The parameter 51o is the 
amplitude value of the atom-field dipole coupling and depends on the position of an atom inside a cavity. To 
treat the nonlinear dynamics of this system Prants et. al. derived first the quantum equations of motion for 
the operators external atomic operators, P and x and slowly varying amplitudes of the field and spin operators: 
a(t) = ae-tof*, a+(t) = a + e iu f\ S ± (t) = S±e^ iut and S z {t) = S z . 

Taking the averages of all operators over an initial quantum state(which is a product of the translational, 
electronic, and the radiation field states) the quantum equations of motions can be replaced by the equations 
for the expectation values of the operators . The equations of motions for these averages are written as 

x = ap 
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p = —2{a x s x + a y s y )sinx 
s x = — Ss y + 2a y s z cosx 

s y = ds x — 2a x s z cos 2 (6) 

Qjjfr — COS IZ/ 



Ciy — &x COS 3j . 

where the expectation values are defined as 

x = kf(x),p= (p)/hk f 
s x = (s_ + s+)/2, s y = (s_ - s + )/2i 

«i = (a + a + )/2, a a = (a - a + )/2i, (7) 
a = Tife? / mi7 Q , r = fl Q t 

In the ref. Q the dynamics of the cavity+atom system at T = is treated by solving Eq. ([6]) and analyzing the 
solutions in terms of Poincare surface sections and Lyapunov exponents. The extension of the zero temperature 
results [9( to non-zero temperature is done here by using the formalism of TFD. 



IV. CAVITY QED AT NON-ZERO TEMPERATURE 

Applying TFD prescription to the cavity QED Hamiltonian (jSJ), we have 

H = ^ + hu f (a+(J3)a(/3) - a+{fi)a{fi)) 

2m 

+m Q [{a(p) sinh0 - a{0) sinh(9 + a+(/3) cosh (9 - a + (/3) cosh(9)S_ 

+ (a+(/3) sinh0 - a+(/3) sinh0a(/3) + cosh<9 - a(fi) cosh 6)S + ) cos(k f x). (8) 

Repeating the same steps as those used in Q for T = we have the temperature-dependent equations of 
motion for the expectation values of coordinate, momentum, spin, annihilation and creation operators 

dx _ P - P 
dt m 
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dP 

— = hk f n [(a + (/3) cosh (9 - a+(P) cosh6> - a[fi) sinh (9 + a(/3) sinh0)S_ 
+(a(j8) cosh 6> + 5+ ((3) sinh - a+ {fi) sinh 9-a{p) cosh 0)S + ] sin fcys 
— ^- = i(u)f - oj a )S+ + 2iCl S z \a(P) sinh 6 - 5(/3)sinh0 + a + (/3) cosh 6* - a(/3) cosh0] cosk f x 

= -i(uj f - uj a )S- - 2iVt S z [a + {(3) sinh0 - a + ((3) sinh + a(/3)cosh0 - a(/3)cosh0] cosfc/i] 



da+ 



-iClo(S- sinh0 — S+ cosh 9) cosfc/i 



= — iOo(<S_ cosh0 — S+ sinh 0) cos kfX (9) 



da 
di 



-i(lo(S+ sinh0 — SL cosh0) cos kfX 



iflo(S+ cosh0 — S- sinh0) coskjx 



where 



dS 

—j^- = {SIq{B\S- — B 2 S + ) coskfx, 



B x = a(/3) sinh 9 - a{(3) sinh + a+(/3) cosh - a + (/3) cosh 0, 



#2 = a + (/3) sinh0 - a + (/3) sinh0 + a(/3) cosh0 - a(/3) cosh<9. (10) 

These equations describe the time evolution of the dynamical variables (expectation values of the operators) 
in the presence of coupling to a thermal bath with constant temperature, T. We note that interaction of the 
spin degrees of freedom with a thermal bath is not taken into consideration in these equations. 



RESULTS AND DISCUSSION 



We have solved numerically the system of equations (j9|) and plotted Poincare surfaces of section (PSS) at 
various values of temperature. In Fig. 1 the Poincare surface section are plotted for: (3 — 2 (a); (3 = 6 (&); 
/3 = 10 (c); 13 = 12 (d). 
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It is clear from these plots that for high temperatures the dynamics is fully chaotic, while by decreasing T 
transition to mixed and regular regime of motion can be observed. To compare the approach for T = and our 
approach for finite T in Fig. 2, PSS are plotted by solving the set of equations (j9)) and (fT0|) (T 0). For T 
about zero our results for PSS agree with those by Prants jsj. 

Following the prescription developed in [9 1 fo r exploring of instabilities in the QED cavity, another character- 
istics of chaoticity, so-called Levy flights [23l - l25| is analyzed. As is well known, chaotic motion in classical system 
have several quantitative and qualitative characteristics, such as phase-space trajectories, Lyapunov exponent 



and Levy flights 24j, 2^ |. The latter are the pieces that appear in the trajectory of a particle in a transition 



from a regular to a chaotic regime of motion. In other words, Levy flights are the chaotic pieces interrupting 
regular behaviour of the trajectory of an oscillating or regularly moving particle 0, Q]. In Fig. 3, Levy flights 
for various values of j3 are plotted. It is clear from this figure that for (a) T = and (b) j3 — 100 the plots are 
similar, with the same number of flights. However, by increasing the temperature (fi = 10 and /3 = 5) leads to 
increasing number of flights. 

To make our treatment more comprehensive we should consider also the behaviour of the maximum Lyapunov 
exponent at different temperatures. The maximum Lyapunov exponent characterizes the mean rate of the 
exponential divergence of initially close trajectories and serves as a quantitative degree of deterministic chaos 
in the system. In Fig. 4 the maximum Lyapunov exponent is plotted as a function of detuning parameter, 
5 at different values of temperature. Again, one can observe "more chaos" in the case of finite temperature 
compared to /3 = 100 case. For higher temperatures (smaller /3) the Lyapunov exponent becomes more higher 
than that for lower temperatures. 

In Fig. 5 the maximum Lyapunov exponent versus atom field detuning 6 and initial atomic momentum po 
is plotted for T — and /3 = 0.5. This plot also shows that increasing of the heat bath temperature leads to 
increasing of maximum Lyapunov exponent for all the values of S and po. However, some "islands" near = 
still exist in this plot. This means that Lyapunov exponent remains as small near S = at T ^ 0. 

In all cases (including the case of T ^ 0) one can observe that A becomes equal to zero at 6 = that means 
becoming of our system integrable for 6 = 0. 

Finally, Fig. 6 presents maximum Lyapunov exponent versus 6 and /3. Again increasing of A for higher T 
(smaller /?) can be observed. Therefore besides the control parameters a and 6, in the case of finite temperature 
we have an additional parameter for controlling the dynamics of the atom in a cavity temperature, T. 



VI. CONCLUSION 

Thus we have studied finite-temperature nonlinear dynamics of an atom coupled to a single mode of the cavity 
field. Applying the formalism of a real-time finite-temperature field theory to the Jaynes-Cummings Hamiltonian 
and using the same approach as that used in [9( we have studied classical dynamics of the " atom+cavity mode" 
system coupling to a thermal bath. The equations of motions for the classical dynamics are obtained by 
" mapping" of quantum dynamics onto classical one as in the ref . [9( . 
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Using the temperature-dependent equations of motion, dependance of the dynamics on heat-bath effects or 
finite temperature effects are considered. The results show that the dynamics is quite sensitive to the small 
changes of temperature. Qualitative characteristics of chaoticity of the system are explored for different values 
of (5. In particular, projections of the Poincare surface section are plotted for different values of temperature. It 
can be seen by comparing the sensitivity of Poincare surface sections, Levy flights and the maximum Lyapunov 
exponent to the changes of /? with those of 5 that the dynamics is more sensitive to the changes of temperature 
than that of S. This implies that the temperature of a thermal bath can be considered as an additional control 
parameter for the dynamics of an atom coupling to a cavity mode. 
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FIG. 1. Projection of the Poincare sections at finite temperature on the plane of the atomic momentum p in units hkf 
and the position in units of kj 1 . (a) P = 2 ; (b) p = 6; (c) /3 = 10; (d) ft = 12. S = 1.92 and s z (0) = -0.863 in all cases. 
x,p are dimensionless. 
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FIG. 2. Projection of the Poincare sections at finite (b) and zero (a) temperature on the plane of the atomic momentum 
p in units hkf and the position in units of kj 1 . (a) T = 0; (b) /3 = 20. 8 = 1.92 and s z (0) = —0.8660254 in both cases. 
x,p are dimensionless. 
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T=6 

£=25 




FIG. 4. (Color online) The maximum Lyapunov exponent A in units of maximal atom-field coupling rate fio versus the 
atom-field detuning S in units flo at zero and finite temperatures: ft = 25; ft = 1; j3 = 0.1; ft = 0.01 and s z (0) = in all 
cases. 
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FIG. 5. The maximum Lyapunov exponent A (in units of maximal atom-field coupling rate Qo) versus the atom-field 
detuning 8 (in units Qq) at different temperatures: T — 0; and (3 — 0.5; (s z (0) = in both cases). 
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FIG. 6. The maximum Lyapunov exponent A (in units of maximal atom-field coupling rate Qo) versus the atom-field 
detuning 8 (in units fio) and /3; (po = 2 and s z (0) = 0). 
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